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A toy model for magnetic connection in black hole (BH) accretion disc is discussed 
based on a poloidal magnetic field generated by a single electric current flowing around 
a Kerr black hole in the equatorial plane. We discuss the effects of the coexistence 
of two kinds of magnetic connection (MC) arising respectively from ([T|) the closed 
field lines connecting the BH horizon with the disc (henceforth MCHD), and ^ the 
closed field lines connecting the plunging region with the disc (henceforth MCPD). 
The magnetic field configuration is constrained by conservation of magnetic flux and a 
criterion of the screw instability of the magnetic field. Two parameters A and a™ are 
introduced to describe our model instead of resolving the complicated MHD equations. 
Compared with MCHD, energy and angular momentum of the plunging particles are 
extracted via MCPD more effectively, provided that the BH spin is not very high. It 
turns out that negative energy can be delivered to the BH by the plunging particles 
without violating the second law of BH thermodynamics, however it cannot be realized 
via MCPD in a stable way. 
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1 INTRODUCTION 

Recently, much attention has been paid to the magnetic con- 
nection (MC) of a rotating black hole (BH) with its sur- 
rounding accretion disc (Blandford f999; Li 2000a; Wang, 
Xiao & Lei 2002, Wang et al. 2003, hereafter W02 and 
W03, respectively) . The MC process can be regarded as one 
of the variants of the Blandford-Znajek (BZ) process pro- 
posed near three decades ago (Blandford & Znajek 1977), 
which involves the closed magnetic field lines connecting a 
BH with its surrounding disc. This mechanism has been used 
to explain a very steep emissivity in the inner region of the 
disc, which is consistent with the XMM-Newton observation 
of the nearby bright Seyfertl galaxy MCG-6-30-15 (Wilms 
2001; Li 2002a; W03). 

Blandford (2002) described several electromagnetic 
ways to extract the rotational energy associated with a spin- 
ning BH, in which the MC between the plunging region and 
the disc (hereafter MCPD) is included. However, MCPD has 
not been discussed in detail in the previous work, and the 
origin of these magnetic field configurations remains unclear. 

Recently, Li (2002b, hereafter L02) discussed the MC 
between the BH horizon and the disc (hereafter MCHD) by 
assuming a toroidal electric current I flowing on a circle of 
r = r' in the equatorial plane of a Kerr BH, i.e.. 



r 



- (A/A)'/' (a/a^)" S{r - r')5(cos6>). 



(1) 



The metric around a Kerr BH is given in Boyer-Lindquist 
coordinates (MacDonald and Thorne 1982, hereafter MT82) 



= - (^1 - 2Ml 

4aA'Ir sin' 



dtdifi + !^dr^ + p^de^ + A^i^dv?^ 



(2) 



where the concerned Kerr metric parameters are defined as 



p-fa') -a'Asin'e 
A = r' a' - 2Mr, vj = (^/p^) 
a = (p'A /A) Lu = 2aMr/A 



p'-^ = r^ + cos' ( 



1/2 



sin 6, 



(3) 



In equation ([Sjl a is the specific angular momentum of the 
BH, which is related to the BH mass M, angular momentum 
J and spin a» by a = J/M = a^M. 

Not long ago, Li (2000b, hereafter LOO) discussed a 
scenario of extracting energy from a Kerr BH through the 
plunging region, where the open magnetic field lines connect 
plasma particles with remote loads. It is argued in LOO that 
the energy extracted from the particles can be so large that 
the particles have negative energy as they fall into the BH, 
if the magnetic field is strong enough. 

Motivated by the above works, we propose a scenario 
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of extracting energy from a Kerr BH based on the mag- 
netic field configuration arising from a toroidal electric cur- 
rent / flowing at the innermost stable circular orbit (ISCO). 
The magnetic field configuration contains both MCHD and 
MCPD, being constrained by the conservation of magnetic 
flux and a criterion of the screw instability of the magnetic 
field. It is shown that MCPD is more important than MCHD 
in transferring energy and angular momentum to the disc for 
low BH spins. It turns out that negative energy can be deliv- 
ered to the BH by the plunging particles without violating 
the second law of BH thermodynamics, however it cannot 
be realized via MCPD in a stable way. 

This paper is organized as follows. In § 2 we present 
a description of our model, in which MCHD and MCPD 
coexist in a large-scale magnetic field generated by a toroidal 
electric current flowing at ISCO. In § 3 we calculate and 
compare the powers and torques in MCHD and MCPD. In 
§ 4 the condition of negative energy of plunging particles 
is studied in based on some reasonable constraints: (i) a 
positive accretion rate, (ii) an increasing BH entropy, (iii) 
a stationary accretion and (iv) a stable energy extraction 
for the plunging particles. In § 5 we compare the efficiencies 
of releasing energy arising from disc accretion, MCHD and 
MCPD. Finally, in § 6, we summarize our main results, and 
discuss some issues related to our model. Throughout this 
paper the geometric units G = c — 1 are used. 



2 DESCRIPTION OF OUR MODEL 

2.1 Magnetic field configuration produced by a 
toroidal current 

It is assumed that the magnetic field in the neighborhood 
of the BH and its surrounding disc is ideally conducting 
and force-free, and the magnetic field, electric field and 
other quantities are measured by zero-angular-momentum 
observers defined by Bardeen, Press & Teukolsky (1972). 
The magnetic field configuration produced by a toroidal cur- 
rent flowing at the inner edge of the disc as shown in Figure 
1, which contains both MCHD and MCPD. In this paper 
we assume that the toroidal current flows at ISCO, being 
represented by the symbol "(8)" . The quantities rH and r^s 
are respectively the radii of the BH horizon and ISCO, and 
they read (Novikov & Thorne 1973) 



= M (l + 



(4) 



r„s = MxLs, 

Xms = {3 + ^2 ± [(3 - Ai) (3 + yli 2A2)f'^f^^ , 

= 1 + (1 - alf' [(1 + a.)^/« + (1 - a,f'] , 
A2 = {ial^AlY'\ 

The quantity re is the radius of the infinite redshift surface 
and it reads 



, {6) = m{i + y/l-al cos^ 6^ 



(6) 



The radius re increases with the polar angle 9, attaining its 
maximum at the equatorial plane, i.e., (7r/2) = 2M. The 
range of the ergosphere is given by 




Figure 1. Magnetic field configuration corresponding to MCHD 
and MCPD, where the symbol "(Xi" represents a toroidal current 
at the inner edge of the disc. 



and the range of the plunging region is 



(8) 



The magnetic flux through a surface bounded by a circle 
with r =const and 6 =const is 

^{r,e-r')^2-KA^{r,e-r'), (9) 

where A^ is the toroidal component of the vector potential 
produced by the toroidal electric current J at r = r' in the 
equatorial plane of the Kerr BH. Defining f = r/M and 
Bo = 2J/M, we have * oc oc Ml = BoM^/2 and 

t2\ 



* (f , 6- ?) = 2-kA^ (r, e- r) / {BqM'' 



(10) 



Equation pup provides a mapping relation between any 
two circles {f\,9\) and (f2,&2) embedded in the same mag- 
netic surface as follows, 

*(fi,ei;r') =*(f2,e2;r'). (11) 

Znajek (1978) and Linet (1979) derived the relation be- 
tween the vector potential A^p (r,6;r') and the current / as 
follows. 



ra sin^ - a sin^ 9 cos 9F2 



-a^sm''9cos9^Fi-rsm'^eF2 + ^0^F'i } 



ra sin^ 9-^F4 - a sin^ 9 cos 9F5 



(12) 



sin^ 9 cos 9 ^Fi-r sin^ gFg + ^6 > 



where in equation p2|l the functions Fi — Fq are defined 
by 

1 =P/(u)P,(cos0), 



r„ <T <T, [9) , 



(7) 



Fi 

F2 
F3 

F4 
F5 



i^P,{u)P; {cos 9), 
1 -,Pl{u)Pl {cos 9).^ 



1 



= Q; {u) Qi (cose), 

{u) Q[ {cos 9) 
=Q[ {u) Q[ (COS0). 



(13) 
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In equation (I13|) we have it = (r — M) / \/ M"^ — a^, 
Pi (z) and Qi {z) are Legendre functions with P/ (2) = 
dPi {z)/dz and Q'l (z) = dQi {z)/dz. The coefficients a[, Qj, 
/Jf and PI are defined as follows: 

(I) for r < r', PI = PI = for all /; but 

«r^ ,.,"!wt.y',,. (^X'\'aP{0)Q[ {u'), (14) 



, ^ 7r(2l + l)I f A'\ ^''^ 

x[- (r'2 + a^)P/(0)Q, (u') 

I r'A' 1 



(15) 



' P/(0)Q;(n')l; 



PI 



(2) for r > r', q[ = ai = for all Z; but 
7r(2/ + l)/ /A'^'''^ 



Z(; + 1)(M2 



^(l^) ^''^^'(0)^/(-'). (16) 



(17) 



ai _ 7r(2i + l)J / 

- i(i + l)(M2-o2) \^ A' y 

xh (r'2 + a')P/(0)Pi(M') 

where in equations (I14|) — p7[) we have A' = A (r = r') and 

A' = A(r = r', 61 = n/2). 



2.2 Mapping relations for MCHD and MCPD 

For the toroidal current located at rms the mapping relation 
between the circle at the BH horizon with the longitudinal 
angle 6 and the circle at the thin disc with the radius r can 
be written as 



* (f , 7r/2; f^,) = * (f^ , 6; fms) ■ 



(18) 



It is assumed that the toroidal current located at rma 
has a very small circular section of radius r'^ = evms, so that 
the poloidal magnetic field is proportional to the radius r' for 
< r' < r^. In this way the infinite magnetic field is avoided 
as rclose to rms , keeping magnetic field configuration outside 
unchanged. In this paper we take e — 10"'* in calculations 
(The influence of the choice of e is discussed in § 6). 

The innermost magnetic surface for MCHD is defined 
as a surface connecting the horizon at 6 = 7r/2 with the disc 
at fo ^ Tms, and it can be determined by 



* (fo, 7I-/2; fms) = * (r^ , vr/2; fms) ■, 



(19) 



where fo = ro/M is the radial coordinate of the disc plane 
intersected with the innermost magnetic surface in MCHD 
as shown by the solid field line in Figure 1. 

Equation (ITT|| can be applied to MCPD by taking Oi — 
O2 = 7r/2. The radial coordinate r^j^ in the plunging region 
and r in the disc are related by 



,7r/2;f„,) = *(f,7r/2;f™s) , 

' PL ^ ^ms, 
fms <f<fo. 



* (f 

r 



PL ' ' 

< fp 



(20) 



Equation (|2Up is the mapping relation for MCPD, which 
relates the radius r-p^ in the plunging region to the radius 
r at the disc. 

The magnetic field configuration can be constrained by 
the screw instability, which will occur if the toroidal mag- 
netic field becomes so strong that the magnetic field line 



turns around itself about once (Kadomtsev 1966; Bateman 
1978). Recently, Wang et al. (2004) argued that the screw 
instability would occur in a magnetized accretion disc, and 
the criterion can be expressed as 



{2-ktx,^/L)B'o/BI < 1, 



(21) 



where and BJ) are respectively the poloidal and toroidal 
components of the magnetic field on the disc, and L is the 
poloidal length of the closed field line connecting the BH 
with the disc. The symbol -cuj^ is the cylindrical radius on 
the disc and it reads 

= ^J~^l Pd 



(22) 



where ^ = r/rms is a dimensionless radial parameter. The 
poloidal component of the magnetic field B^ can be written 
as 

P [Ap/Apf^^ d1>(r,V2;r^,) 

^ 2^ Tr • (23) 

Equation (|23p is derived based on the calculation on the 
magnetic fiux between two adjacent magnetic surfaces, i.e., 



d*(r,7r/2;r„ 



~Bo 2Try/grrgipipdr. 



(24) 



In W02 the powers and torques in the BZ and MCHD 
processes are derived based on an equivalent circuit consist- 
ing of a series of loops, which correspond to a series of adja- 
cent magnetic surfaces. By using Ampere's law the toroidal 
magnetic field B'f, can be expressed as 



B 



2/f,g/(. 



(25) 



The electric current Imc flowing in each loop of the equiv- 
alent circuit for MCHD is given by 

HD _ Ae^ + Ae^ 
- AZH + AZn- 

The resistance AZd in equation (|26|l can be neglected 
due to the perfect conductivity of the disc plasmas. The 
quantities Ae^ and Ae^ are electromotive forces due to 
the rotation of the BH and the disc, respectively, and they 
read 

Ae„ = [(A*)^/27r] Qh, 
Ae„ = - [(A*)^/27r] Qo, 

where (A^/)^ and (A*!*)^ are the fluxes sandwiched between 
two adjacent magnetic surfaces threading the BH horizon 
and the disc, respectively. The minus sign in the expres- 
sion of As^ arises from the direction of the fiux. We have 
(A^*)^ = (A*!")^ and the two can be calculated by using 
equation (|24p . The quantities Qh and Q,d are respectively 
the angular velocity of the BH horizon and that of the disc, 
and they read 

^ (28) 



(27) 



2Mf, 



fin = 



M(f3/2 4- a.) 

Incorporating equations (^1])- 
rion of the screw instability as follows 



we have the crite- 



{Ad/AdY^'' {-dZn/dr) {uj^a] 



9 = ir/2 



< L/{2nuj^) 



(29) 



2{nH -no) 

In equation (|29p dZn is the resistance of the BH horizon 
between the two adjacent magnetic surfaces, and we have 
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dZn/dr = {dZn/de) {dO/dr) 



(30) 



where d9/dr < can be calculated by using the mapping 
relation (fTSt . 

Incorporating equations (I2ip — (|30|) . we can derived the 
critical radius of the screw instability, fg = r^/M. As 
was demonstrated by Uzdensky (2005), the pressure of the 
toroidal magnetic field wound up by the rotating BH can- 
not be compensated in the outer regions of the magneto- 
sphere. Thus we regard the radius as the outer boundary 
of MCHD region as shown in Figure 1. The critical radius 
is related to the angle 9s by the mapping relation (|18p . 
More specifically, the mapping relation for MCHD can be 
expressed as 



Os <e< 7t/2, 

ro < r < r „ . 



(31) 



Thus the magnetic field configuration produced by the 
toroidal current is constrained by the conservation of mag- 
netic flux and the criterion of the screw instability. By using 
equations @, ®, (|5J, (fT^ and (^5)) we have the curves of 
the radii , fms, fe, fo and versus a, as shown in Figure 
2. 



It is shown that the radii r. 



fo and r„ decrease 



monotonically with the increasing a,, while the radius fe 
remaining a constant "2" in the equatorial plane. It is easy 
to check from equations ^ and ^ that fms = re holds for 
a» = ^8/9 = 0.9428, beyond which the inner edge of the 
disc will enter into the ergosphere. 



3 POWERS AND TORQUES IN MCHD AND 
MCPD 

The motion of the plunging particles is very difi'erent from 
that in the disc. The velocity of the radial inflow in the 
disc is much smaller than the Keplerian velocity, while the 
radial velocity of plunging particles is comparable to the Ke- 
plerian velocity. The motion of plunging particles becomes 
even more complicated in the presence of large-scale mag- 
netic fields. In this paper, instead of resolving MHD equa- 
tions in the plunging region, we introduce a parameter A 
to describe the MCPD effect on the motion of the plunging 
particles. 



3.1 Angular velocity of plunging particles 

It is reasonable to assume that the accreting matter mainly 
consists of electrons and protons. Considering the toroidal 
motion of the accreting electrons and protons in the plung- 
ing region and the magnetic field configuration in Figure 1, 
we infer that the electrons of negative charge and the pro- 
tons of positive charge are accelerated inwards and deceler- 
ated outwards by Lorentz force, respectively. Thus a current 
is driven outwards in the plunging region by Lorentz force. 
This situation is analogous to those in the BZ process and 
in MCHD, where a current is driven in the longitudinal di- 
rection at the horizon due to the BH rotation relative to a 
poloidal large-scale magnetic field (MT82; W02; W03). And 
we expect to derive the power and torque in MCPD by using 
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Figure 2. Curves of f^, fms, fe, fg and versus a* in dot- 
dashed, dashed, thick solid, dotted and thin solid lines, respec- 
tively, for < a, < 0.998 (a) and 0.8 < a, < 0.998 (b). 



an equivalent circuit analogous to that given for MCHD in 
W02 and W03. 

Starting from the Lagrangian of a test particle (Shapiro 
& Teukolsky 1983), we can derive the expression for the 
angular velocity Slpz, of the particle in the plunging region 
as follows. 



npL = 



(f^^ -2) {Lms/M) + 2a,Er, 



M {f %^ + alr^^ + 2ai) E,ns - 2a, {Lms/M) ' 



(32) 



where Ems and Lms are the specific energy and angular 
momentum of the accreting particles at ISCO, and they read 



Ejjis — 



1 2Vms ^ a,Tms 



1 3r^s -4" 2(1, 



1/2 ' 



Lms/M = 



( 

' ms \ 



1 - 2a,fJJ'' + a*f;js 



1 St'ttt.s -|- '2(l:^V')jis 



1/2 



(33) 



(34) 



Considering the magnetic extraction of energy and an- 
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gular momentum, we suggest that the angular velocity, fip^ , 
of the plunging particle is related to Q.pl by 



f^PL = (1 — A cost) VlpL, 



(35) 



where A is a parameter to adjust the ratio of fi™^ to Q.pL, 
of which the value range is < A < 1. The parameter r is 
defined as 



-f„)^/2 



(36) 



Incorporating equations (|32|) . p5[) and H36p . we find 
that the angular velocity JlpL is decelerated at most 
to {l — X)Q,pL at fpj^ — {fma+fjj)/2, and it connects 
smoothly with Q.h and Q.d a,t and rma, respectively. 

For the given values of A the angular velocity fip^ is a 
function of a« and fp^^ , and we have the curves of fipj, /Qo 
varying with f^j^ for the given BH spins as shown in Figure 
3. 

It is shown in Figure 3 that the profile of f^p^ is af- 
fected significantly by the parameter A at the middle of the 
plunging region, where a less fipj^ corresponds to a larger A. 
The value of A is crucial for MCPD in transferring energy 
and angular momentum, and we shall discuss its value range 
based on some physical considerations. 



3.2 Derivation of powers and torques in MCHD 
AND MCPD 

We can derive the expressions for power and torque in 
MCPD by using an analogous equivalent circuit given in 
W02, and the electromotive force due to the rotation of the 
plunging matter in the plunging region is written as 



(37) 



where (A'I')p^ is the magnetic fiux threading the plunging 
region between the two adjacent magnetic surfaces and it 
reads 



(38) 



The electric current in each loop of the equivalent circuit 
can be written as 



PD AepL + Ae 

J^MC — 



where Aej-, and AspL are expressed by equations (|27|l and 
(|37p . respectively. The power and torque transferred from 



the pluuj 


jing region 


A T>P° 










AZpL ) 


^J-MC — 


A 



\2tv J 



(40) 

PL \ ■2n ) ^ 



(41) 



AZpL 

In equation (|39[) the resistance of the disc plasma is 
neglected due to the perfect conductivity. The resistance of 
the matter in the plunging region, AZpL , can be determined 
by the requirement of continuity of the electric current in 
the two adjacent loops bounded by the innermost magnetic 
surface for MCHD, i.e.. 



(42) 



Incorporating equations ((23, ([33), (EH and 

we have the ratio of AZpl to AZh as follows, 



AZf 



AepL+Asr 



AZe 



Ae, 



Ae, 



(43) 



Considering that the innermost magnetic surface for 
MCHD is located at Vp^ = r^, we infer that Qpl — > i^H 
due to the frame dragging of the Kerr BH. From equation 
(|43|) we have A^p_t ~ AZh with Aepi, — Ae^ for the two 
equal adjacent magnetic fiuxes at ro. 

As a simple analysis we assume that the unknown sur- 
face resistivity of the plasma fiuid is equal to that of the BH 
horizon, i.e.. 



4tv. 



Thus the resistance AZpL in equation 
as 

/^drp. 



dZ, 



PL 



R 



(44) 

can be written 



(45) 



PL ' PL -F 



-2f, 



Finally, the total power and torque in MCHD are ob- 

MC '-^^MC 



tained by integrating AP"fn and AT^rn over the MCHD 



region as follows. 

P5§/Po = JT^ [BE if)] ' {MQh - MQd) (M^n) 
xf{e,f;at) {-df/de) df, 

tZc/To = [BE (f)] ' {MQh - MQu) 
xf{e,f-a:,) {-df/de) df, 



(46) 



(47) 



where the function / {9, f; a«) in equations (|46[) and (|47[) are 
give by 



/( 



20 



^<l^„2 -2 



(^2 csc2 e + ^l-aj-l^ (f2 + aj-2i=) sin 9 



(48) 



Similarly, the total power and torque in MCPD can be 
-f^gandAT^.g 

region as follows. 



calculated by integrating AP[jc and ATfjc over the MCPD 



P[w/Po = IZ^ [BE (r-)]' (Affl?J^ - Amo) {M^d) 
^g{fpL^f\o-*) {''df/dfp^)df, 

Tmc/To = J'l [BE (r)] ' {Mn^^L - MQo) 
^g{fpL^f\o-*) {~df/dfp^)df, 



(49) 



(50) 



where the function g {rpj^,f;a,) in equations (|49p and (|50[) 
are give by 



gifp 



^{1 + aifpl + 2fplal) {f%^ +ai- 2fp,) (51) 



f'^+a'i—2f 



In equations (|46|) — (|50|) we have —df/dO > 0, 
-df/dfp^ >0, Po = BlM"^, To = B^M^ and 

BE if) . BE /Bo = - (^-y-)^'^ '^^'^ -/^^ . (52) 



dr 



6 




Figure 3. The curves of flpi^/Qo and $7p^/Qo versus with a» = 0,0.3594,0.9428 in (a), (b), and (c), respectively, where Clpi^/Qo 
is plotted in solid line with A=0, and Q^j^/Qq is plotted in dashed and dotted lines with A=0.1 and 0.3, respectively. The value range 
of fp^ is fjj < fp^ < fms and SIq = l/M = 2.03 X 10^^ (r ad/ s)(M/Mq)~'^ 
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(b) 

Figure 4. (a) The curves of -P^g/-Po (solid Une) and -Pjgg/Po 
(dashed line) versus a*; (b) the curves of T^j^/Tg (solid line) 
and Tjgg/To (dashed line) versus a». The parameter A = 0.5 is 
assumed. 



By using equations (ge)), (gT]), (ggjl and dSOj we have the 
curves of P^g / Po, T^g/To, Pmc/Po and rf.g/To versus 
a, with a given value of A as shown in Figure 4. 

Inspecting Figure 4, we find the following results: 

U]) For the increasing a, with the given A, -Pjv/c 
1'mc decrease monotonically in a very fast way, while Pmc 
and Tmc increase monotonically in a much slower way. 

@ The quantities Pmc and T[jc are greater than P"jg 
and Tmc for low BH spins, and the former two become 
negative for high BH spins, while Pmq and T^g remain 
positive for all BH spins. 

Incorporating equations (|49p — (|51[). we have the con- 
tours of Pf,g = and P[i^/P^ig = 1, T^^g = and 
TMc/'^mg = 1, in a, — A parameter space as shown in 
Figures 5 and 6, respectively. 

From Figures 5 and 6 we find that the BH spin a, and 
the parameter A are two crucial parameters for MCPD, and 
it is more effective than MCHD in transferring energy and 
angular momentum into the disc for a* and A in the region 
III. On the other hand, MCPD could be less effective than 
MCHD in this aspect for the two parameters falling into the 
very narrow region II as shown in Figures 6b and 7b. The in- 
tegrated MCPD transfers of energy and angular momentum 
are reversed for the two parameters in region III. 



4 BH EVOLUTION AND NEGATIVE ENERGY 
IN PLUNGING REGION 

4.1 BH Evolution equations 

Based on the conservation of energy and angular momentum 
we have the evolution equations of a Kerr BH by considering 
disc accretion with MCHD and MCPD as follows, 

dM/dt = EmsMo - Pug - Pmc, (53) 

dJ/dt = LmsMo — T^g — Tmc- (54) 

Incorporating equations (|53p and H54p . we have the evo- 
lution equation for the BH spin as follows, 

da,/dt = (L„,Mo - T^,g - T[,^) 

-2M-'a. {E^sMn ~ Pftg - PlFc) ■ ^ ' 
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Figure 5. The contours of Pj^g = (solid line), P[,c:/Pmc = ^ 
(dashed line) in a, — A parameter space for (a) < a, < 0.9428 
and (b) 0.835 < a, < 0.880. We have P^g < < Pj(Jg, < 
Pmc < ^MC ^'^'^ < P/fg < PjPg in regions I, II and III, 
respectively. 
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Figure 6. The contours of Tj^g = (solid line), T^^/^MC = ^ 
(dashed line) in a* — A parameter space for (a) < a* < 0.9428 
and (b) 0.845 < a, < 0.890. We have T|^g < < Tjgg, < 
^MC < '^MC ^'^'l < Tj»g < TjPg in regions I, II and III, 
respectively. 



Incorporating equations (|53[) — (|55p . 
change rate of the BH entropy as follows, 



we derive the 



THdSn/dt = dM /dt - flndj/dt = {E„,s - 0.hL^s) 



PD\ 
MC 



(56) 



where Th and Sh are the temperature and entropy of the 
Kerr BH, and they read (Thorne, Price & MacDonald 1986) 



Th 



47rAf(l + g) ' 



Sh ^2TvM^{l + q). 



(57) 



Considering the energy and angular momentum trans- 
ferred between the plunging region and the disc, we think 
that disc accretion should be affected significantly by 
MCPD. There must be some relations between the large- 
scale magnetic field and the accretion rate Md- As a mat- 
ter of fact these relations might be rather complicated, and 
would be very different in different situations. One of them 
is given by considering the balance between the pressure of 
the magnetic field on the horizon and the ram pressure of 



the innermost parts of an accretion flow (Moderski, Sikora 
& Lasota 1997), i.e., 



{B^y/ (Stt) = Pram ~ pc"" ~ Md / {47Vrl) . 



(58) 



Considering the fact that the poloidal magnetic field 
arises from a toroidal current at ISCO, we replace and 
in equation (|58p by the magnetic field near ISCO and 
, respectively. In addition, we introduce a parameter am 
to adjust the ratio of the magnetic pressure to the ram pres- 
sure of the accretion flow in a thin disc, and assume the 
relation between the accretion rate and the poloidal mag- 
netic field at ISCO as follows. 



Mn = a„ 



[B^nirms)frl,s/2+{MD) 



MC ■ 



(59) 



where the second term at RHS of equation (|59p is the MCPD 
correction to the accretion rate and it reads 



(Md) 



(60) 



Substituting equations (|28|) and (|41|) into equation (|60p . we 
have 



2M[Bg (f )]" (f (ng^ -na)9(f .f ;a,) (df /df ) 



(61) 



Equations (|60|l and (|6ip are derived based on the con- 
servation of angular momentum. It is noted that the sign 
of df/dfp^ is negative, and the MCPD correction to the 
accretion rate depends on both the sign and magnitude of 

Since (Md)^^^ depends on the disc radius, the accre- 
tion rate expressed by equation (|59p varies with the disc 
radius, which might give rise to an unstable accretion disc. 
Fortunately, the viscous force can adjust itself to compensate 
for the excessive transfer of angular momentum by magnetic 
torque and keep the accretion rate constant. How to con- 
strain the concerning parameters to ensure a stable accre- 
tion disc? We shall deal with this issue based on the transfer 
of angular momentum in the disc. 

The flux of angular momentum H is transferred from 
the plunging region to the disc, and it is related to T^/^ by 
A-nrH = OT^'jq j dr. Based on the conservation of energy and 
angular momentum Page & Thorne (1974) derived the rate 
of angular momentum transferred by the viscous torque in 
the disc as follows. 



Qvis — 



i-dflD/dr) 



-fMl 



(62) 



where and are the specific energy and angular mo- 
mentum of the accreting particles, respectively. The quantity 
M^'^ can be regarded as the first term at RHS of equation 
((59)) . and the function / = ArKrF j M"^" is related to the 
radiation flux F given by Page & Thorne (1974). 

Being required by a stationary accretion, the rate of 
angular momentum transferred by the viscous torque should 
be no less than that due to MCPD, i.e.. 



gvia ^ 47rrff = QtI^ I dr. 



(63) 



Based on equation (|62|l and (|63|) we have the condition 
for a stationary accretion as follows. 



Facc > 1, 

FA.c=g.rsl{dTl°cldr), 



(64) 



where the function Facc is evaluated at f = (fms + ) /2- 
Another constraint to the accretion rate is given based 
on a stable extraction of energy. The accretion rate must 
be great enough to keep a stable magnetic extraction of en- 
ergy from the plunging particles. As a rough estimate, the 
decreasing rate of the plunging particles' kinetic energy can 
be written as 



(65) 



which should be no less than the MCPD power at any po- 
sition of the plunging region. Considering the fact that the 
plunging particles' energy cannot be extracted completely 
via MCPD, we suggest the condition for a stable energy ex- 
traction as follows, 



where Pj(/c ('^pi, ) the MCPD power integrated to the po- 
sition at fp^ — {rms +ru)/2. Thus the condition (|66|) can 
be rewritten as 



Fsne 



4.2 



Pk 



PMc{r^^)^0. 



(67) 



Condition for Negative Energy of Plunging 
Particles 



Not long ago, LOO discussed the possibility of negative en- 
ergy carried by accreting particles onto a BH due to the 
magnetic extraction of energy to the remote loads. An in- 
teresting issue is whether negative energy can be realized 
via MCPD. We are going to discuss this possibility with the 
constraints based on several physical considerations. 

Following LOO, we have the condition for the negative 
energy of the plunging particles as follows. 



PpL = MnErr 



(68) 



Pk ^ PlFc 



(66) 



where Ppl is the rate of the net energy brought into the BH 
by the plunging particles. The negative energy is determined 
by PpL < with the following constraints. 

Il]) The rate of the BH entropy given by equation (|56p 
is not negative, which implies that the BH entropy can never 
decrease required by the second law of BH thermodynamics 
(Thorne, Price & MacDonald 1986); 

([2]) The modifled accretion rate expressed by equation 
(|59|) remains positive to provide the plunging particles in a 
continuous way; 

(131) The stationary accretion is guaranteed by equation 

dMl); 

Q A stable energy extraction from the plunging parti- 
cles is ensured based on equation (|57)) . 

Incorporating equations (|56p — H68[l. we have the con- 
tours of dSn/dt = 0, Md = 0, Facc = 1, Fsnergy = aud 
PpL = in the a, — A parameter space as shown in Figure 
7. 

In Figure 7 the shaded region indicated "Negative En- 
ergy" is bounded by four contours, i.e. Ppl ^ 0, dSn/dt ^ 
0, Md ^ and Facc ^ 1 in dashed, dotted, solid dot-dashed 
lines, respectively. However, the negative energy is excluded 
by the contour of Fsnergy ^ in thick solid line, above 
which the stable energy extraction from the plunging par- 
ticles is required. The region above the thick solid line and 
the dot-dashed line is indicated "Stable Accretion and En- 
ergy Extraction" , in which energy extraction via MCPD can 
work in a stable way. 

Thus we conclude, although the negative energy can 
be delivered to the BH via MCPD with the validity of the 
second law of the BH thermodynamics, it cannot be realized 
in a stable way in our model. 



5 EFFECTS OF MCPD AND MCHD ON 
EFFICIENCY OF ENERGY RELEASE 

Although MCHD and MCPD are related to the magnetic 
field configuration supported by a toroidal current flowing on 
the equatorial plane of the Kerr BH, these two mechanisms 
are different in several aspects. 

(Ul Energy and angular momentum are extracted and 
transferred from the Kerr BH to the disc by virtue of MCHD, 
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Figure 7. The contours of PpL = (dashed line), dSu/dt = 
(dotted line), M^i = (solid line), F/^^^ = 1 (dot-dashed line) 
and FEne-rgy = (thick solid line) in a, — A parameter space 
with < a* < 0.9428 and < A < 1 for (a) = 0.039 and (b) 
Q,„ = 0.030. 
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spectively in solid, dashed and dotted lines for (a) am = 0.039, 
A = 0.55 and (b) a„ = 0.03, A = 0.6. 



while those are extracted and transferred from the plunging 
particles to the inner disc by virtue of MCPD. 

((2]) MCHD can work in a suspended accretion state as 
suggested by van Putten & Ostriker (2001), while MCPD 
cannot operate without continuous replenishment of accret- 
ing particles. 

Therefore a stable accretion and energy extraction are 
indispensable in the model containing the coexistence of 
MCHD and MCPD. Based on the BH evolution equation 
(|53p . we can easily obtain the efficiencies of releasing energy 
as follows. 

r? = 1 - {dM/dt) /Md ^VoA+ Vmc + Vmc , (69) 

where rj^A, Vmc s-nd Vaw ^re respectively the efficiencies 
due to disc accretion, MCHD and MCPD, converting the 
rest energy of the accreting particles into the radiation en- 
ergy. These efficiencies are expressed as 



{VoA = 1 — Ems, 
VMC = PMg/Mn, (70) 

Incorporating equation (|70p with the concerning expres- 
sions for Ems, Mn, Pmc ^-^i^ Pmc, have the curves of 
VdA, tIaic 8-nd t/mc versus a, for the different values of q™ 
and A as shown in Figure 8. 

The efficiency v^/c is greater (less) than rjfjQ for the 
lower (higher) BH spin, while it is generally less than rj^A 
as shown in Figure 8. In order to compare the efficiencies 
rj^/c ^-iid r/j-jA in an accurate way we have the contours of 
Vaic /VdA = 1, Facc = 1 and FEnergy = iu fl, — A param- 
eter space as shown in Figure 9. 

Inspecting Figure 9, we find that efficiency jyff^ is gen- 
erally less than 77^ a in a wide value range of a* and A for 
the given values of am, and the possibility of > Vd^ 

is confined in a very narrow shaded region indicated by 
''VMc/vnA > 1" in Figure 9. 
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Figure 9. Contours in a» — A parameter space of V^i^/VoA = 
li F^Energy = and -F^cc = 1 in dotted, solid and dot-dashed 
lines, respectively. The parameter am=0.039 and 0.03 are given 
in Figures 10a and 10b, respectively. 



6 DISCUSSION 

In this toy model the magnetic field configuration is gen- 
erated by a single electric current fiowing at ISCO of a 
thin disc, and we regard this as a starting point for dealing 
with the magnetic field configuration generated by the elec- 
tric current distributed in more realistic way around a BH. 
We concentrate the discussion on the magnetic extraction 
of energy and angular momentum from the plunging parti- 
cles, and we describe our model phenomenally by using two 
parameters instead of resolving the extremely complicated 
MHD equations within ISCO. The parameter A is used to 
describe the angular velocity of the plunging particles, and 
the parameter am is invoked to adjust the accretion rate af- 
fected by MCPD. It is noted that the parameters a,, A and 
Om are crucial for a stable MCPD process, although the 
MCPD power and torque are independent of the parameter 
am according to equations (P^ - (f5T]) . 

The key point in this model is how to constrain these 
parameters based on some reasonable considerations. In this 
model four requirements are given: (i) a positive accretion 



rate, (ii) an increasing BH entropy, (iii) a stationary ac- 
cretion and (iv) a stable energy extraction for the plunging 
particles. Required by these constraints the following results 
are obtained in our model. 

(m MCPD could be more important than MCHD in 
transferring energy and angular momentum to the disc. As 
shown in Figures 4-6, the power and torque in MCPD are 
stronger than those in MCHD for low BH spins. 

(121) Although negative energy can be delivered to the BH 
by the plunging particles with the increasing BH entropy it 
cannot be realized via MCPD in a stable way as shown in 
Figure 7. 

The efficiency of releasing energy due to MCPD could 
be greater than that due to MCHD for low BH spins, while 
it is generally less than the efBciency of disc accretion for a 
wide range of a* and A as shown in Figure 9. 

It is assumed in § 2.2 that the toroidal current lo- 
cated at rms has a very small circular section of radius 
= ETms, and e = 10~* is taken in calculations. The pa- 
rameter e is introduced to avoid an infinite magnetic field 
as rclose to the toroidal current. The magnetic field out- 
side ri not affected by the value of £, while the 
contribution of the magnetic field to the MCPD effects can 
be neglected if e is small enough. Thus the influences of e 
on the MCPD effects are very little. For example, required 
by Pmc/Pmc = Iwith a* =0.294, we have A=0.6187569089, 
0.6187569092 and 0.6187569342 for 10"^ 10"" and 10~^ 
respectively. The variation of A arising from that of e is 
AA/Ae < 3 X 10"^ 

Another issue related to our model lies in prescription 
(1441) . in which the surface resistivity of the plunging region 
is assumed to be Rpl = Rh = 47r. The resistivity Rh is 
defined such that the electric field at the BH horizon is equal 
in magnitude to the magnetic field. We can show that the 
prescription (|44p does not lead to the result that the electric 
field at the surface of the plunging region is greater than the 
magnetic field there. 

The toroidal component of the magnetic field at the 
surface of the disc is expressed by equation (|25p . Similarly, 
we can derive the expressions for the toroidal magnetic field 
at the surface of the plunging region as follows, 

bIl = (^J7PLa)e^^/2- (71) 

The toroidal electric field is zero in a stationary axisym- 
metric magnetosphere (MT82), while the poloidal electric 
field at the surface of the plunging region is related to 
by Ohm's law, i.e., 

aEpL = RplIIFc/2hzupl, (72) 

where aEpL is the electric field converted to a 'per unit 
globe time' basis (MT82). Incorporating equations (|25p and 
^72)) we have 

Epl/b'^l = RpL/'iTv. (73) 

Thus the prescription (|44p leads to EpL ~ Bpj^ and EpL < 

BPL = V^B^J^+CB^J". 

In this paper we do not fit the very steep emissivity 
observed in some BH systems, such as galaxy MCG-6-30- 
15. In fact, a very steep emissivity can be also produced by 
virtue of MCPD, since the energy and angular momentum 
transferred from the plunging region are very concentrated 
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in the inner disc. We shall work out a more realistic model 
to relate MCPD to the observations in our future work. 
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